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We discuss the problem of trapping and mobilization of nonwetting fluids during immiscible 
two phase displacement processes in porous media. Capillary desaturation curves give residual 
J~^ ' saturations as a function of capillary number. Interpreting capillary numbers as the ratio of viscous 

to capillary forces the breakpoint in experimental curves contradicts the theoretically predicted force 
balance. We show that replotting the data against a novel macroscopic capillary number resolves 
the problem for discontinuous mode displacement. 
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PACS: 47.55Mh, 81.05Rm, 61.43Gt, 83.10Lk, (47.55Kf) 



I. INTRODUCTION 



A displacement of one fluid by another within a porous medium poses challenging problems of micro-to-macro-scale 
transitions which have received considerable attention from physicists in recent years |l|-|2(]|]. For reviews the reader 
is referred to p^ , pl[ |. Apart from the omnipresent quenched correlated disorder and structural heterogeneity, fluid- 
fluid and fluid-solid interactions generate metastability and hysteresis phenomena which have resisted quantitative 
prediction and understanding. A central problem of great practical importance is the prediction of residual nonwetting 
fluid saturation after flooding the pore space with immiscible wetting fluid |g^-Q. Obvious interest in this problem 
arises from enhanced oil recovery |25y26[| or in situ remediation of soil contaminants P7| , p8[ . 

Microscopically the laws of hydrodynamics governing the pore scale processes are well known. The complexity of 
the microscopic fluid movements and the lack of knowledge about the microstructure and wetting properties, however, 
renders a detailed microscopic treatment impossible. Instead one has to resort to a more macroscopic treatment. 

Desaturation experiments p9|-p6| show that the conventional macroscopic description is incomplete gS^M. We 



shall begin our discussion by reminding the readers of this incompleteness. We then exhibit another problem |37|. This 
arises from the fact that microscopic capillary numbers seemingly cannot represent the force balance in a desaturation 
experiment. 

Given these problems our main objective is to show that the recent analysis of |^ gives the correct force bal- 
ance between macroscopic viscous and capillary forces for continuous mode desaturation experiments. For so called 
discontinuous mode displacements it leads to bounds on the size of residual blobs. 



II. PROBLEMS OF TWO-PHASE FLOW EQUATIONS 



Let us begin our discussion with the standard macroscopic equations of motion for two phase immiscible displace- 
ment. Macroscopic equations of motion describe multiphase flow on length scales large compared to a typical pore 
diameter. Hence they are applied to laboratory samples with linear dimensions on the order of centimeteres as well as 
to whole reservoirs measuring kilometers and more. Consider a system having lengths L^:, Ly, L^ in the three spatial 
directions. The pore space is assumed to be filled with two immiscible fluids denoted generically as water (index w) 
and oil (index o). The equations read [^-|o|J3^ 

</>^ = V.(^^[VP„-pw5 0e]| (I) 

^dS^ - ^^^'°[V(Pw+Pc)-Po5 0e]l (2) 



dt y iio 

and they are supplemented with the constitutive relationships 

kl,{^,t) = ¥y,{Sn{^,t)) (3) 

fc;(x,i) = fc;(^„(x,t)) (4) 

P,{i^,t) = F,(5„(x,t)) (5) 

where 

g„(x,t)^ f-("/^'^- . (6) 



The variables in these equations are the pressure field of water denoted as Pw, and the water saturation S^. The 
saturation is defined as the ratio of water volume to pore space volume. Pressures and saturations are averages 
over a macroscopic region much larger than the pore size, but much smaller than the system size. Their arguments 
are the macroscopic space and time variables (x, i). The saturations obey S^i < Sw < 1 — Sor where the two 
numbers < Swi,Sor < 1 are two parameters representing the irreducible water saturation, SvM^^ and the residual 
oil saturation, Sor- The residual oil saturation gives the amount of oil remaining in a porous medium after water 
injection. The normalized saturation Sn varies between and 1 as Sw varies between Swt and Sor- The permeability 
of the porous medium is given by the absolute (single phase flow) permeability tensor K. The porosity 4> is the volume 
fraction of pore space. The two fluids are characterized by their viscosities /^vv, A*o and their densities Pw,Po- The 
terms pg O e represent the gravitational body force where e^ = (0,0,-1) is a unit row vector pointing along the 
negative z-axis, and g is the acceleration of gravity. The orthogonal matrix 

cos Uy sin ay 
O = ( 1 I (7) 

— sin ay cos ay 

describes an inclination of the system. Here a rotation around the y-axis with tilt or dip angle ay was assumed. 
The macroscopic capillary pressure Pc is defined as the pressure difference between the oil and the water phase. 
The constitutive relation for Pc assumes that the capillary pressure function Pc depends only on the saturation 
pl[ . In addition to Pc the dimensionless relative permeabilities are assumed to be functions of saturation only, 
^w ~ K^ASw ^k'^^ = ^ k'^{Sw). They represent the reduction in permeability for one phase due to the presence of the 
other phase ||3^,|39|. The three constitutive relations fc^(5'w), fco(5'w) and Pc{Sw) are assumed to be known from 
experiment. Equations (El) and (0) arc coupled nonlinear partial differential equations which must be complemented 
with large scale boundary conditions. For laboratory experiments the boundary conditions are typically given by a 
surface source on one side of the sample, a surface sink on the opposite face, and impermeable walls on the other 
faces. For a geosystem the boundary conditions depend upon the well configuration and the geological modeling of 
the reservoir environment. 

The problem with the macroscopic equations of motion (|l|)-(^ arises from the experimental observation that the 
par ameters Sor and Swi are not constant and known, but depend strongly on the flow conditions in the experiment 
|p9|- p6| , p6[ . Hence they may vary in space and time. More precisely, the residual oil saturation depends strongly on the 
microscopic capillary number Ca = Pwv/crow where w is a typical flow velocity and aow is the surface tension between 
the two fluids. The capillary desaturation curve ^^^(Ca) is shown in Figure for unconsolidated glass beads and 
sandstones [^^. Such curves contradict clearly to the assumption that the functions kl^{Sw),ko{Sw) and Pc{Sw) 
depend only on saturation. Instead they show that kl^{S„),ko{S„) and Pc{Sw) depend also on velocity |4|,Q and 
pressure. Hence they depend on the solution and cannot be considered to be constitutive relations characterizing the 
system. The dependence shows that the system of equations of motion is incomplete. 

The breakpoint in a capillary desaturation curve marks the point where the viscous forces, which attempt to 
mobilize the oil, become stronger than the capillary forces, which try to keep the oil in place. For this reason the 
capillary desaturation curves are usually plotted against Ca which represents the microscopic force balance between 
viscous and capillary forces. From Figure |l| it is seen that the theoretical force balance corresponding to Ca = 1 on 
the abscissa and the experimental force balances represented by the various breakpoints at Ca <C 1 differ by several 
orders of magnitude. Plotting residual saturation against the correct ratio of viscous to capillary forces should result 
in a breakpoint at Ca « 1 [^^'E^ ■ 

We now proceed to show that a partial understanding of the macroscopic force balance can already be obtained 
from the traditional equations of motion. The main result of this analysis is a preliminary experimental validation of 
the dimensional considerations introduced in |37ll. 



III. DISCONTINUOUS VS. CONTINUOUS MODE DISPLACEMENT 



Before embarking on a discussion of the force balance during capillary desaturation it is crucial to emphasize an 
important difference in the way capillary desaturation curves are measured. 

1. In the first method of measuring S'or (Ca) the oil in a fully oil saturated sample is displaced with water at a 
very low Ca. After the oil flow at the sample outlet stops and several pore volumes of water have been injected 
without producing more oil the flow rate is increased. Again the oil flow is monitored until no more oil appears 
at the outlet. In this way the flow rate is increased iteratively until Sor has fallen to zero. After the first 
injection the oil configuration is discontinuous (or disconnected) and it remains disconnected throughout the 
rest of the experiment. This mobilization mode will be called discontinuous. 

2. In the second method of measuring capillary desaturation curves one again starts from a fully saturated sample, 
and performs a waterflood at a given value of Ca. After oil flow ceases at the outlet the residual saturation is 
determined. Then the sample is again saturated fully with oil, a new value of Ca is chosen, and the injection 
is repeated. In this experiment the injection starts always with a connected oil phase contrary to the previous 
experiment where the oil phase is discontinuous at higher Ca. This mobilization mode will be called continuous. 

The capillary number required to reach a given Sor is known to be much lower in the continuous mode than in the 
discontinuous mode JS^JS^I . This is also seen in Figure |l| which shows continuous mode data and discontinuous mode 
data for the unconsolidated glass beads. From the fact that the equations of motion are only valid in the subinterval 
S„i < S„ < 1 — Sar it follows that they cannot be applied to discontinuous mobilization mode experiments. They 
should however be valid for the continuous mode to the extent that the equations themselves are correct. Therefore 
we discuss next the macroscopic force balance predicted by these equations for the continuous mobilization mode. 



IV. MACROSCOPIC FORCE BALANCE 




It was shown in |37| that the balance of macroscopic viscous and capillary forces is not represented appropriately 
by the microscopic capillary number Ca, and that a new macroscopic capillary number Ca should be used instead. 
Here we generalize those results to anisotropic and inclined porous media and apply them to replot the Sor data 
obtained for the continuous mode displacement. To this end we rewrite the equations (|l|) and (y) in dimensionless 
form. Introducing the matrix 



(8) 

and defining I and L through 

I = (detL)(i/3) (9) 

L = / L (10) 

one defines dimensionless quantities 

x = LX = /Lx (11) 



V = L-^V = (/L)-iV. (12) 

Define k and K through 



i^xx i^xy i^xz 

iV ^ I kxy rCyy Kyz | \^'^) 

i^xz f^yz i^zz 

fc=(detK)(i/3) (14) 

K = fcK (15) 

The time scale is normalized as 

t^hihMl (16) 

using the volumetric flow rate Q. Time is measured in units of injected pore volumes. The pressure is normalized 
using the macroscopic equilibrium capillary pressure as [}44| , p7[ 

P^PbP (17) 

where 

Pb = Pc aSy,, - Sor + l)/2) (18) 

is the pressure at an intermediate saturation. PdSyv) is the equilibrium capillary pressure function. 
With these definitions the dimensionless macroscopic capillary numbers for oil and water, defined as 

give an expression of the balance between macroscopic viscous and capillary forces. The macroscopic gravity numbers 

Pwy "■' 
&, = &„i^^ (22) 

Mw Po 



express the viscous to gravity force balance. Finally the gravillary numbers 



CjIw — 



Ph 



(23) 



LtLo Cjiv 



Po 



(24) 



express the ratio between gravitational and capillary forces. The well known bond number, measuring the magnitude 
of buoyancy forces, is given as 



Bo = Glw - Gl, 



W ^J^O 



(25) 



in terms of the gravillary numbers. 

With these definitions the dimensionless equations of motion may be rewritten as 



dt 



-{ 



Afc,': 



^^w ^ Pw ^ Gr^^ g 



} 



(26) 



dS„ 
dt 



= v -lAk:^ 



[a kl [Ca^'v(Pw + Pc) - GrJ g] } 



(27) 



where the dimensionless matrix 



A = L^KL^ 



~Epr L^Lyk L^L^k 

I'^h I'^h l-^k 

I r\jxy I' "^yy ^ yz 

LxLyk Lfk LyL^k 






l^k 



yz 



I- kzz 



(28) 



\ LxLzk LyLzk Ljk / 



contains generalized "aspect ratios" . The vector 



(LOe) 



L, 



■ sinQ!j^, 0, 



L, 



■ cos a. 



(29) 



represents the effect of dip angle and geometric shape of the system on the gravitational driving force. 



V. APPLICATION TO EXPERIMENT 



We are now in a position to plot the capillary desaturation curves against the macroscopic capillary number Ca 
which represents the balance between macroscopic viscous and capillary forces. To do so we have searched the 
literature for capillary desaturation measurements and found Refs. [p2 25 29|-pl|,E5|-E8|,p5l . Unfortunately none of 



the publications contains all the necessary flow and medium parameters to calculate Ca. Measuring all the flow 



parameters for a displacement process is costly and time consuming, and hence they are rarely available (see also 
POj). While permeability, porosity and the fluid parameters such as vicosities and surface tensions are usually 
available capillary pressure data, relative permeabilities and residual saturations are not routinely measured. Hence 
we extract the required parameters from different publications assuming that they have been measured correctly and 
are reproducible anywhere and at all times. In spite of all the uncertainties it is well known that the capillar y p ressure 
curves of unconsolidated sands with various grain sizes can be collapsed using the Leverett-j-correlation ||50|,^,y . 
The Leverett-j-correlation states essentially that 



Pc{Sw) = (^owy-rJiSw)- (30) 

Often the formula contains in addition an average contact angle at a three phase contact. We do not include the 
wetting angle as it is generally unknown, and including it would not change our results significantly. Similar to 
the capillary pressure data, the Sor data for unconsolidated sands seem to be well established [p4|36| in spite of 



larger fluctuations of the results. Because most Pc- and 5or-data are available for unconsolidated sand and standard 
sandstones such as Berea or Fontainebleau we limit our analysis to these two cases. 

The experimental S^r data analyzed here are taken from Ref. [ p2[ for sandstones and from Rcf. p& for uncon- 
solidated glass beads. In Figure |l| we show the capillary desaturation curve for oil- water displacement in a typical 
sandstone (sample No. 799 from ||3^) using star symbols. These data were obtained in the continuous mode of 
displacement. The values of the surface tension and fluid viscosities for this experiment are given in Table |. We also 
show continuous mode S'or-data for unconsolidated glass beads from Figure 6 in |Q • To calculate Ca we have used 
the values given in Table |. 

Plotting the data against Ca we obtain Figure 0. It is seen that the continuous mode displacements give a breakpoint 
for Ca w 1 while the discontinuous mode displacements have their breakpoint at a higher value. This is consistent 
with the idea that the traditional equations of motion (m) and (0) should be applicable to the continuous mode but 
not to the discontinuous case. 

The result obtained here is consistent with the theoretical predictions from p7| . To fully validate our use of Ca as 
a correlating group for plotting continuous mode Sor data, however, it would be desirable to vary Ca by varying the 
system size /. If the Sor curves obtained for different media and length scales I also show their breakpoint at Ca « 1 
this would give further evidence for the applicability of the traditional equations of motion. We consider it possible, 
however, that deviations appear indicating a breakdown of the equations also for continuous mode displacement 
11,11 ■ 

As stated above the equations of motion are not applicable to discontinuous mode displacements because of the 
constraint Swi < Sw < 1 — Sor- Nevertheless we can use the group Ca to estimate an upper bound for the size 
of residual blobs. As the flow rate is increased the residual blobs whose size is so large that the viscous drag forces 
on them exeed the capillary retention forces will break up and coalesce with other blobs which may again break up 
and coalesce further downstrea m ||3^ , |33[ | . The condition that the viscous forces dominate the capillary forces Ca > 1 
predicts that after a flood with Ca the porous medium contains only blobs of linear size smaller than 

Iblob < „ 1 ■ (31) 

Pbk 

This result is of importance for microscopic models [ p2| of breakup and coalescence during immiscible displacement. 
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TABLE I. The sample and flow parameters used in the calculations. The permeabilities for the glass bead packs were 
obtained using the relation between their radius and permeabilities published in |p5[ on page 47. The capillary pressure was 
estimated using the Leverett-j'-function correlation (pd). For sandstone the j-function from Ref. [pll was used and for bead 
packs the one from Ref. [pOl was used. The j-function was evaluated close to {Swi — Sor + l)/2 according to (|l8|). 
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FIGURE CAPTIONS 

Figure 1: Experimentally measured capillary desaturation curves (capillary number correlations) for bead packs 
[|36| (solid lines with circles and squares) and sandstone ||32| (star symbols) as a function of microscopic 
capillary number Ca — fJ-m^/'^ow The dashed line is the desaturation curve for continuous mode 
displacement. The dash-dotted line marks the plateau value for sandstone. 

Figure 2: Same as Figure n^ but plotted against the macroscopic capillary number Ca — iiwQ/{Pbkl) from cq. 
(Q9f). Note that the breakpoint for continuous mode displacement occurs around Ca « 1. 
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